Composite boson many-body theory for Prenkel excitons 



Monique Combescot (1) and Walter Pogosov (1,2) 
(1) Institut des NanoSciences de Paris, Universite Pierre et Marie Curie, CNRS, 140 rue de Lourmel, 75015 Paris and 

(2) Institute for Theoretical and Applied Electrodynamics, 
Russian Academy of Sciences, Izhorskaya 13/19, 125412 Moscow 
(Dated: January 30, 2009) 

We present a many-body theory for Frenkel excitons which takes into account their composite 
nature exactly. Our approach is based on four commutators similar to the ones we previously 
proposed for Wannier excitons. They allow us to calculate any physical quantity dealing with N 
excitons in terms of "Pauli scatterings" for carrier exchange in the absence of carrier interaction and 
"interaction scatterings" for carrier interaction in the absence of carrier exchange. We show that 
Frenkel excitons have a novel "transfer assisted exchange scattering", specific to these excitons. It 
comes from indirect Coulomb processes between localized atomic states. These indirect processes, 
commonly called "electron-hole exchange" in the case of Wannier excitons and most often neglected, 
are crucial for Frenkel excitons, as they are the only ones responsible for the excitation transfer. We 
also show that in spite of the fact that Frenkel excitons are made of electrons and holes on the same 
CO ■ atomic site, so that we could naively see them as elementary particles, they definitely are composite 

objects, their composite nature appearing through various properties, not always easy to guess. The 
present many-body theory for Frenkel excitons is thus going to appear as highly valuable to securely 
tackle their many-body physics, as in the case of nonlinear optical effects in organic semiconductors. 
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PACS numbers: 71.38.-y, 71.35.Aa 



I. INTRODUCTION 



Even if nature makes things somewhat more complex, we can roughly say that there are essentially two types of 
excitons: Wannier excitons [1| and Frenkel excitons Q]. 
\ (i) Wannier excitons are found in conventional (inorganic) semiconductors. They are constructed on valence and 
<— ] ■ conduction extended states. Their Bohr radius, controlled by Coulomb interaction between electrons and holes, 
screened by the semiconductor dielectric constant, is as large as a few tens of interatomic distances. This makes their 
O . binding energy rather small. 

(ii) Frenkel excitons are the excitons of organic semiconductors. They are constructed on highly localized atomic 
, states. Their "Bohr radius", controlled by the atomic wave functions, is of the order of a single molecular block. This 
^ ' makes their binding energy as large as a fraction of e.V. or more. 
I/"") , Wannier excitons made of two fermions, the conduction electron and valence hole, were up to now commonly 
considered as elementary bosons, the possible fermion exchanges being included through effective scatterings, by lack 
of appropriate procedure. Over the last few years, we have developed this missing many-body theory [1, H[, through a 
procedure which allows us to treat the composite nature of Wannier excitons exactly. This theory now has its specific 
diagrammatic representation in terms of "Shiva diagrams" [j| which not only allow to visualize the tricky iV-body 
exchanges which exist in these systems but also to calculate them readily. This new many-body theory has been 
successfully applied to a variety of physical effects such as Faraday rotation @ and oscillation Q , spin precession 8] 
and teleportation 0, polariton and exciton Bose-Einstein condensation [Icj . In each of these effects, the composite 
nature of the excitons plays a key role; in particular this new theory allows us to show that the replacement of 
Wannier excitons by elementary bosons misses the correct detuning behavior of all nonlinear optical effects induced 
by nonabsorbed photons. 

Since organic materials are of potential importance for new electronic devices, the proper handling of many-body 
effects with Frenkel excitons along a theory similar to the one we have developed for Wannier excitons, is of high 
current interest. For readers not aware of this new composite exciton many-body theory, its fundamental ideas are 
briefly explained in a short review paper ||. More details can be found in an extended Physics Report [1] and the 
references it contains to more than 30 published works. 

In order to tackle Frenkel excitons on a firm basis we have, in a previous paper [llj . rederived the basic concept 
leading to these excitons, starting from first principles, i.e., a microscopic Hamiltonian written in first quantization, 
with electron and hole kinetic energies, electron-ion, electron-electron, and ion-ion Coulomb interactions (This ion- 
ion term, which does not affect the electron motion, makes the system at hand neutral, in this way insuring the 
convergence of the overall Coulomb contribution in the large sample limit). Through a grouping of terms appropriate 
to highly localized electron states - this grouping being different from the one leading to valence and conduction bands 
on which the Wannier excitons are constructed - we make appear the atomic states of each ion. When overlaps between 
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relevant atomic states of different ion sites are negligible, these states can be used as a one-body basis to rewrite the 
system Hamiltonian in second quantization. In a second step, we transform the ground and excited atomic creation 
operators into electron and hole creation operators. Among the various Coulomb terms which keep the number of 
electron- hole pairs constant, one is insuring the local neutrality of the excitations. The resulting set of lowest excited 
states, made of these locally neutral excitations, with one electron and one hole on the same site, are degenerate in 
the absence of Coulomb interaction between sites. We then introduce these intersite interactions. Among them, one 
is a transfer term which destroys one electron-hole pair on one site and recreates it on another site. This transfer 
term splits the degenerate subspace made of one-site excitations into linear combinations of these excitations which 
are nothing but the Frenkel exciton states. 

Frenkel excitons, made of electron-hole pairs localized on all possible atomic sites, are composite bosons, definitely. 
Their many-body effects thus have to be handled with extreme care if we want to fully trust the obtained results. 
This is why it can appear as appropriate to construct for them a full-proof procedure similar to the one we have 
used for Wannier excitons [!, [3, [l^, HU [H| > these Wannier excitons also being linear combinations of electrons and 
holes but in extended states - not in localized states as for Frenkel excitons. Actually, due to this high localization, 
we can be tempted to believe that carrier exchanges between Frenkel excitons should reduce to zero; so that these 
excitons should behave as elementary bosons. In the present work, we show that this idea is actually a very naive 
one: Frenkel excitons do behave as elementary bosons for some quantities, but differ from them in other important 
respects, making their many-body physics not easy to guess. This is why the many-body theory we here propose is 
going to be highly valuable because it allows us to tackle many-body effects with Frenkel excitons in a fully secure 
way. 

Before going further, let us point out one important difference between Wannier and Frenkel excitons: Wannier 
excitons are controlled by one length for both, the interaction and the extension of the relative motion wave function. 
This length is the Bohr radius induced by the screened Coulomb interaction between electron and hole and calculated 
with the effective masses of the conduction and valence bands. The situation for Frenkel excitons is far more complex. 
We could at first think that Frenkel excitons have two lengths: the atomic wave function extension for the electron-hole 
pair "relative motion" and the Bohr radius for Coulomb interaction. Actually, Frenkel excitons have not one but two 
relevant Coulomb interactions: one for direct Coulomb processes between sites; the other for indirect processes, this 
last interaction being the one responsible for the excitation transfer between sites. In addition, as the electrons and 
holes on which the Frenkel excitons are constructed are not free, the proper effective mass to use in a "Bohr radius" 
is a priori unclear. Finally, the tight-binding approximation underlying Frenkel exciton corresponds to reduce the 
atomic wave function extension to zero, so that the length associated to the electron- hole "relative motion" de facto 
disappears. Due to these intrinsic difficulties and the expected interplay between these various "physical lengths" , the 
many-body physics of Frenkel excitons is going to be far more subtle than the one of Wannier excitons, controlled by 
one length only. 

The goal of the present paper is to settle the basis for a composite boson many-body theory appropriate to these 
Frenkel excitons. Using the formalism presented in the present work, we are going, in the near future, to address to the 
calculation of various observables, such as the ground state energy of N Frenkel excitons, the scattering rate of such 
excitons, and the nonlinear susceptibility of materials with Frenkel excitons, the later problem being very important 
for applications. Possible extension of the present theory can be to include effects of interaction with the vibrational 
modes, e.g., disorder and exciton-phonon coupling. This can be done by adding the corresponding potential to the 
Frenkel exciton Hamiltonian. 

The paper is organized as follows. 

In Section II, we recall the main steps leading to Frenkel excitons, as derived in our first work on these excitons 

In Section III, we derive the "Pauli scatterings" for carrier exchange between two Frenkel excitons, starting from 
the " deviation-from-boson operator" of these excitons [HI, [l3j] . 

In Section IV, we introduce Shiva diagrams for carrier exchanges and show their similarities with Shiva diagrams 
for Wannier excitons [j| • We also point out an intrinsic difficulty of these Frenkel exciton Shiva diagrams linked to 
the tight-binding approximation which makes atomic electron-hole states used to represent Frenkel excitons, not a 
clean complete basis for all electron-hole pair states. 

In Section V, we use these Pauli scatterings to calculate some relevant scalar products, including the normalization 
factor N\Fn for N identical Frenkel excitons. We find that, as for Wannier excitons [H|,[l5|], Fn which plays a key role 
in all exciton many-body effects, decreases exponentially with N due to carrier exchanges between composite bosons. 
This Fn turns out to have an extremely simple compact form due to the fact that Frenkel excitons are characterized 
by one quantum number only, the center-of-mass momentum Q, while Wannier excitons in addition have a relative 
motion index v. This is going to greatly simplify all results on many-body effects linked to the composite nature of 
these particles. We also discuss the physical meaning of the relevant dimensionless parameter rj which controls carrier 
exchanges between excitons. We show that, in spite of our first understanding in the case of Wannier excitons, this 
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key parameter for many-body effects is not linked to the spatial extension of the electron-hole relative motion wave 
function but to the total number of cxcitons the sample can accommodate - otherwise rj would reduce to zero in the 
case of Frenkel excitons. This new physical understanding of rj is one of the important results of the present work; it 
provides a significant step toward mastering the many-body physics of composite bosons, in general. 

In Section VI, we derive the closure relation which exists for iV exciton states and we show that, while the prefactor 
for Wannier excitons [IB] is (1/Nl) 2 , the prefactor for Frenkel excitons is 1/N\ as for N elementary bosons. The fact 
that the prefactor for closure relation of composite bosons made of 2 fermions is not always (1/iV!) 2 , as we found for 
Wannier excitons, is also one of the important results of the paper. More on this prefactor can be found in Ref. [l7l |. 

In Section VII, we calculate the creation potential of one Frenkel exciton, as necessary to properly describe the 
Coulomb interaction of this exciton with the rest of the system [l4j]. This Frenkel exciton creation potential 
actually splits into three conceptually different terms: One term comes from direct and exchange Coulomb processes 
inside the same site responsible for local neutrality; the second term originates from the indirect Coulomb processes 
between sites responsible for the excitation transfer; the third term comes from direct Coulomb processes between 
sites, this last contribution being similar to the one found for Wannier excitons. 

In Section VIII, we calculate the interaction scatterings associated to the three parts of the creation potential. We 
find that the one associated to direct Coulomb processes between sites depends on the momentum transfer between 
the "in" and "out" excitons, as in the case of Wannier excitons [HKlIj]. The scattering associated to the neutrality 
term is physically very strange as it does not depend on the momenta of the " in" and " out" excitons involved in the 
scattering. It however disappears from the physics of these excitons, when summation over all possible processes is 
performed. The third scattering, associated to the excitation transfer between sites, can also be seen as strange as 
it does not depend on the momentum transfer but just on the momenta of the "out" excitons. This is physically 
linked to the fact that the excitation transfer does not induce any momentum change. This third scattering actually 
corresponds to a "transfer assisted exchange" : the coupling between two excitons comes from carrier exchange through 
a Pauli scattering while the energy-like contribution to this scattering comes from the transfer part of the Frenkel 
exciton energy resulting from indirect Coulomb processes. It can be of interest to note that this transfer assisted 
exchange scattering, somehow unexpected at first, has similarity with the "photon assisted exchange scattering" 
between two polaritons we recently identified [Tq |. The symmetry between "out" and "in" states is actually restored 
when considering energy conserving process, the scatterings being then equal, as shown in Section IX. In this Section 
VIII, we also calculate the Coulomb exchange scatterings when Coulomb interaction is followed by carrier exchange. 
We show that, as a consequence of the tight-binding approximation, the Coulomb exchange scatterings associated to 
direct and transfer processes reduce to zero. 

To enlighten the interplay between the Pauli scattering and the various interaction scatterings between two Frenkel 
excitons, Section IX is devoted as an exercise to the calculation of the matrix element of the Frenkel exciton Hamilto- 
nian in the two-exciton subspace. We make clear the part coming from non-interacting elementary bosons, the part 
coming from non-interacting composite excitons and the part coming from interaction between these excitons. We 
also compare this matrix element with the one obtained by assuming that Frenkel excitons are interacting elemen- 
tary bosons. We again show that, in order to have equality between matrix elements for composite and bosonized 
Frenkel excitons, we are forced to use physically unacceptable effective scatterings, i.e., scatterings which would induce 
non-hermiticity in the Hamiltonian. 

In Section X, we discuss the various results obtained on Frenkel excitons by comparing them with similar results 
obtained for elementary bosons and for Wannier excitons. We show that Frenkel excitons have similarities with both. 
This in particular means that they definitely differ from elementary bosons by various means which makes all intuitive 
handling of these excitons rather dangerous. 

Section XI is a brief state-of-the-art. We discuss the approaches proposed by the two main groups who have tackled 
interaction between Frenkel excitons [2(| Hi], H2, HH, 0, [2f| [26[ and point out some of the advantages of the present 
theory. 

In the last Section, we conclude with a brief summary of the main results derived in the present work. 

This paper is definitely quite long: We have chosen to include all fundamental results on the composite boson 
many-body theory for Frenkel excitons in a unique paper, for the reader to have at hand all the necessary tools to 
securely tackle any physical effect dealing with these excitons. Results on the Pauli and interaction scatterings (l3l.[Tl|. 
the normalization factor 

[HI ], the closure relation [l6[ and the Shiva diagram expansion Q were the purposes of 
different papers at the time we were constructing a similar composite boson many-body theory for Wannier excitons. 
Through all these works, we have however learned what should be done to be fast and efficient; This is why we can 
now propose a rather compact presentation for all these important pieces of the composite Frenkel exciton puzzle. 
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II. ALL THE WAY TO FRENKEL EXCITONS 
A. Electron and hole creation operators 

Let us call a* = a\ =1 n the operator which creates one electron in the atomic excited state v = 1 on site n with n 
running from 1 to the total number N s of sites in the periodic lattice. In the same way, the hole creation operator 
&}j = a„ = o. n destroys one electron in the atomic ground state v = on site n. In the tight-binding limit in which 
the products of atomic wave functions are such that <p*,(r — R„/)</?„(r — R„) ~ for all n' ^ n when (v, v') = (0, 1), 
these operators obey 

a n ,a\ + a\a n > = {a„ s 4} - S n 'n - {K'^D (2.1) 

{a n >,a n } = {b n >,b n } = {a„/,^} = ~ {a n >,b n } (2.2) 

They can thus be used as a one-body basis for electrons and holes in problems in which only the ground and first 
excited atomic states play a role. 

B. Electron- hole Hamiltonian 

The Hamiltonian for these free electrons and holes reads 

Heh = e e ^ a li a n + £h ^2 b n h n (2.3) 

n n 

where the electron and hole energies e e and Sh are slightly different from the excited and ground state atomic energies 
due to Coulomb interactions with the "jellium" made of one electron in the atomic ground state of each site n. A 
similar difference also exists for Wannier excitons, the electron and hole energies being slightly different from the 
conduction and valence electron energies. 

The one-electron- hole-pair eigenstates of H e h are a^ b^ \v) (see Fig. 1(a)), their energy being e e + Sh- Since n e 
and n h can run from 1 to N s , these states form a N% degenerate subspace. 

C. Pair Hamiltonian 

If we now take into account the intrasite direct and exchange contribution of the Coulomb potential, the resulting 
pair Hamiltonian is given by 

Hpair — H e h -)- Vintra 

V intra = -£ ^ a^6],6„a„ = -£ B\ B n (2.4) 

n n 

ajj6^ creates one electron-hole pair on site n. The energy eo gained by such a pair reads as 

eo - V R = {I I) - V R = (? I) (2.5) 

The Coulomb couplings between atomic states Vr. (^f correspond to atom on site n\, going from v\ to u[, while 

atom on site 712 goes from v 2 to v' 2 , the two sites being at R„ = R m — R„ 2 . It is related to the ground and first 
excited atomic wave functions <£„(r) through 

V R (:j 2) = J dv.dv^l, (nK, (r 2 ) _ ^ - R| Vv2 (r 2 )^ (n) (2.6) 

Note that Eq. (2.6) makes the two terms of e re& l while the orthogonality of atomic states makes it positive; so that 
Vintra forces local neutrality in the electron-hole excitation - as physically expected due to the electrostatic energy 
cost to separate electron from hole. 
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Consequently, among the N% one-pair states a^jfi \v), the N s states with the electron and hole on the same site 
have the same lowest energy (see Fig. 1 (b)). This leads to 

{Hpair — Epair) \n) = (2.7) 

In) = 46+ \v) = Bt \ v) (2.8) 
where E pa i r = e e + £h — so- Note that, since 4 — 0, due to Eq. (2.2), the electron-hole operators B n are such that 

Bl = (2.9) 

D. Free exciton Hamiltonian 

The next step is to introduce the "electron-hole exchange" part of the Coulomb interaction, i.e., the so-called 
"transfer term". It describes intersite indirect Coulomb processes between v — and v = 1 states and allows 
excitation transfer between sites (see Fig. 1(c)). This leads to the "free exciton Hamiltonian" defined as 

E^X ~ Elpair ~t~ Vtransf 
Vtransf = E Vn ni -R„ 2 (? o) 4 1 4 1 & ™2 a "2 

™1#™2 

= E Vn ni -n n2 {il)BlB n2 (2.10) 

nutn 2 

The prefactor of this transfer term corresponds to Coulomb process in which one electron-hole pair on site ni recom- 
bines, i.e., the atom on site ni is deexcited from v = 1 to v = 0, while the atom on site m is excited from v = to 
v = 1. This transfer term splits the N s degenerate subspace made of states \n) into N s Frenkel exciton states defined 
as 



1 N 3 

I^Q>=^b) = - 71 =E elQ - R "H ( 2 - n ) 



n=l 



While \n) corresponds to a localized excitation, Frenkel excitons \Xq) correspond to delocalized excitations through 
the sum over the R„ sites it contains. By noting that 



Vtransf \Xq) = E E ^-R™ (? I) ^^44 l«> 



= -4= E e iQ - R -44 1«> E e ^ Q - (R -- Rm v Rn _ Rro (s j) (2.i2) 

while the sum over m is independent on n due to translational invariance, we readily see that the \Xq) states are 
such that 

(H ( p - E Q ) \X Q ) = (2.13) 
where the exciton energy Eq is related to the transfer coupling through 

E Q =e e +e h -e +V Q (2.14) 

v Q = E e ~ iQR M?J)=VQ ( 2 - 15 ) 
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Note that Vq is real as expected for an energy since, due to Eq. (2.6), we do have 

V£ (? I) = Vn (J ?) = V- n (? J) (2.16) 

Also note that there is a plus sign in the exponential of the Frenkel exciton creation operator while there is a minus 
sign in the exponential of the Frenkel exciton energy. 

While Wannier excitons are characterized by a center of mass momentum Q and a relative motion index, Frenkel 
excitons are characterized by Q only. This is going to make the many-body physics of Frenkel excitons far simpler 
than the one of Wannier excitons, the "in" excitons (Qi, Q2) and the "out" excitons (Q[, Q 2 ) of a scattering process 
being simply linked by Qi + Q^ = Q1 + Q2, due to momentum conservation, without any additional degree of freedom. 

By noting that the creation of an exciton Q and the creation of an excitation on site n are linked by 

V ly s n =l 

Bt=4&t = 1 £ e -<Q.R.j^ (2.18) 

V s Q 

while due to the lattice periodicity, we do have 

^e iQ '( R n'- R ") = N.6 n , n (2.19) 
Q 

£y(Q'- Q ).R„ =Ns§Q , Q (2.20) 

n 

it is easy to show that the neutrality term Vi ntra in Eq. (2.4) reads as ~£oJ2qBqBq while the transfer term 
Vtrans reads as J2q Vq^q^q. Consequently, the free exciton Hamiltonian can be written as 

=e e J2 + e h J2 blb n + S x (2.21) 

n n 

S x = J2t Q Bl l B Q (2.22) 
Q 

where we have set 

Cq = V Q - e (2.23) 

Note that, while exciton operators can be used to rewrite the Coulomb part Sx of this free exciton Hamiltonian, this 
is not possible for the kinetic energy parts which still reads in terms of fermion operators. Also note that, although 
coming from two-body Coulomb interaction, Sx appears as a diagonal operator in the Frenkel exciton subspace. 
This is going to have important consequence for the scattering of two Frenkel excitons associated to this part of the 
Hamiltonian. 

E. Interacting exciton Hamiltonian 

Coulomb interaction also contains intersite direct contributions coming from electron-electron and hole-hole 
Coulomb interactions, as well as terms coming from direct electron-hole Coulomb processes, the exchange electron-hole 
term between sites being already included in through V trans f. Consequently, the interacting exciton Hamiltonian 
ultimately reads as 

H x = + V coul (2.24) 
Vcoui = V ee + V hh + V e { t ] (2-25) 
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where the three terms of the Coulomb interaction between electrons and holes, shown in Fig. 1 (d, e, f), correspond 
to direct processes. They precisely read 

V ee = l J2 V R^-Rn 2 (I 1) alai 2 a n2 a ni (2.26) 
V hh = \ E VR^-Rnz (o o) blblb n2 b ni (2.27) 

Kf r) = E *V-*„ 2 (J 5) ^W 2 a^ ni (2-28) 

Note that none of these three potentials can be exactly written in terms of exciton operators £>q. Also note that 
the potential V cou i plays a role on states having more than one electron-hole pair since it contains operators acting 
on different sites, while each site can be occupied by one electron-hole pair only - in the absence of spin degrees of 
freedom. The introduction of these spin degrees of freedom and the quite subtle polarization effects they induce, are 
out of the scope of the present paper. 



III. CARRIER EXCHANGES 
A. Free pair operators 

In the previous section, we made appear the excitations of one electron-hole pair on site n through the creation 
operators B\ = a^b^ (see Eq. (2.8)). From Eqs. (2.1, 2), it is easy to show that the commutators of these free pair 
operators are such that 

B n 'B n — B n B n > = [B n ',B n ] = (3.1) 

[B n ,,Bl]_=5 n , n -D n , n (3.2) 
where the deviation-from-boson operator D n > n for free pairs is given by 

D n < n = S n ' n (ala n + blfin) (3.3) 
Its commutator with another pair thus reduces to 

[D n , n ,Bl] =25 n , n 5 nm Bl (3.4) 

This quite simple result physically comes from the fact that electrons and holes have highly localized atomic wave 
functions; so that electron-hole pairs can feel each other by the Pauli exclusion principle when they are on the same 
site only. 

B. Pauli scatterings for Frenkel excitons 

Let us now turn to delocalized excitations, i.e., Frenkel excitons with creation operators given in terms of free pair 
operators by Eq. (2.17). From the commutators of these free pairs, we readily find 

[B Q ,,B Q ]=0 (3.5) 



Bq> , B Q 



= <5q'q - £>Q'Q (3.6) 



where the deviation-from-boson operator for Frenkel excitons is given by 

v (e) , a CO 

*Q-Q' ^Q-Q' 



^Q'Q = Aq-Q' + Aq-Q' (3.7) 
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in which we have set 



A< e > = 



1 iP.R, t 



(3.8) 



and similarly for Ap with a^a n replaced by b^bn. 

Using the expression of free pair operators B\ in terms of Frenkel exciton operators given in Eq. (2.18) and the 
system periodicity through Eqs. (2.19, 20), it is easy to show that 



Dq'q, -Bp 



= — B 



N ^P+Q-Q' 



(3.9) 



If we now compare this equation with the standard definition of the Pauli scattering, deduced from the one for Wannier 
excitons, namely 



^q,4]_ = E{ a P^)+ a (q : p)} b - 



(3.10) 



we are led to identify the Pauli scattering of two Frenkel excitons with 



A 



Q' Q\ _ 
P' p 



= A 



P' Q 
Q' P 



Q', P+Q 



(3.11) 



When compared to the Pauli scatterings of Wannier excitons [lH [l4j]> the one for Frenkel excitons appears as 
structureless: It barely contains the momentum conservation expected for any scattering. This extremely simple form 
of Pauli scattering can be traced back to the lack of degrees of freedom in the electron-hole pairs on which Frenkel 
excitons are constructed. Indeed, electrons and holes are in the two lowest atomic states by construction, i.e., the 
states for which the tight-binding approximation is valid. These \v = (0, 1), n) states do not form a complete set for 
all electron-hole pairs but can be used as a basis for the Frenkel exciton relevant states. In contrast, Wannier excitons 
are constructed on all the extended electron-hole states k e and k/j which have the same center-of-mass momentum 
k e + kft, = Q. In addition to the center-of-mass momentum Q, Wannier excitons thus have a momentum p for the 
relative motion of the electron and the hole in the exciton. This momentum p is related to the electron and hole 
momenta on which Wannier exciton is made, through k e = p+a e Q and k^ = p+a^Q where a e = 1 — ah — 
m e /(m e + m h ). 



C. Many-body effects coming from carrier exchanges 



As for Wannier excitons, in order to derive many-body effects induced by carrier exchanges between any number 
of Frenkel excitons in an easy way, it is necessary to construct commutators similar to the ones of Eqs. (3.6, 9), but 
for N excitons 

moi. They read 



D P , P ,B. 



Q 



! N 

1 n: 



-B, 



Q+P-P'^P 



,B 



t-iV-l 



(3.12) 



tJV 
Q 



NB^- 1 (5 P , Q - D P , Q ) 



N(N-l) ]N _ 2 t 

N D Q i3 2Q-P' 



(3.13) 



as easy to recover by iteration. These commutators are just the ones for Wannier excitons (see Eq. (5) in Ref. |15| ) 
with the Pauli scattering replaced by its value for Frenkel excitons, namely N~ 1 Sq>+p> ,q+p- 



IV. SHIVA DIAGRAMS 



Shiva diagrams for carrier exchanges between Wannier excitons [5] have been shown to be extremely powerful. 
They not only allow to visualize all carrier exchanges which take place between excitons, but also to calculate the 
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physical effects these exchanges induce readily. This makes them as valuable for the many-body physics of composite 
bosons as the Feynman diagrams for the many-body physics of elementary quantum particles. This is why it can be 
of interest to settle similar Shiva diagrams for carrier exchanges between Frenkel excitons. However, as shown below, 
their intuitive handling turns out to be less obvious due to the fact that the B\ \v) states on which Frenkel excitons 
are constructed, do not form a clean complete basis for all one-pair states - in contrast to the free conduction electron 
and valence hole states on which the Wannier excitons are constructed. 

A. Pauli scattering for two Frenkel excitons 

Although possibly surprising at first, the very simple Pauli scattering of Frenkel excitons found in Eq. (3.11) still 
has the same formal structure as the one for Wannier excitons fl3l Ibij ] . Indeed, by writing the wave function of one 
electron on site n as (r \ \v) — y^=i(r — R„) and the wave function of one hole on site n as (r| fet \ v) — ^* =0 (r-R n ), 
where (p u =o,i are the atomic wave functions, the exciton Q wave function reads as 

(r h ,r e |Q) = $q(e» = -j= J] e lQ - R "<^(r e - R„)^(r h - R„) (4.1) 

* 3 n 

where we have set tp e = ip v= \ and ip^ = <p* =0 . If we now calculate the diagram corresponding to hole exchange 
between two excitons shown in Fig. 2(a), it reads in terms of the Frenkel exciton wave functions given in Eq. (4.1) as 

J d{r}$;,(e 1 ,/i 2 )$Q/(e2,/ii)$Q(e 2 ,/i 2 )$p(ei ) /i 1 ) (4.2) 
This quantity contains integrals like 

/^ :( r e -R„0^(r e -R„)-^ (4.3) 

for highly localized atomic states; so that, due to Eq. (2.20), the value of this fermion exchange diagram (4.2) reduces 
to 

_L £ e *(-P'-Q' + P+Q).R„ = _L, p;+Q , p+Q (4 . 4) 

s n 

which is nothing but the value of the Pauli scattering A ^p, p^j obtained in Eq. (3.11). Consequently, as for Wannier 

excitons, the Pauli scattering of two Frenkel excitons can be represented by the Shiva diagram for fermion exchange 
between two excitons, shown in Fig. 2(a). The unique (relevant) difference comes from the fact that, due to the 
lack of relative motion index, the Pauli scattering for Frenkel excitons is the same for an electron exchange or a hole 
exchange. 

B. Carrier exchange between more than two excitons 

If we now turn to the carrier exchange between three Frenkel excitons which corresponds to the Shiva diagram of 
Fig. 2(b) and we calculate it along the standard rules for Shiva diagrams [f|, i.e., by just writing what we read on 
the diagram, we find using Eqs. (4.1, 3) and Eq. (2.20) that 



Q' 3 


Q 3 \ r 


Q 2 


Q 2 =/ 


VQi 


Qi/ J 



$Qi(ei, /ii)$q 2 (e 2 , /i 2 )$Q 3 (e3, h 3 ) 



_ Jp 5 Qi+Q'2+Q'3> Q1+Q2+Q3 
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This result agrees with the decomposition of three-body exchange in terms of a succession of two-body exchanges 
shown in Fig. 2(b), namely 

It is easy to check that this is also true for the four-body exchange shown in Fig. 2(c). And so on... So that the 
Shiva diagram for carrier exchange between ./V excitons just corresponds to 

n n-i S Q'i+---+Q' n > Qi+-+Qjv 

C. Difficulties with Frenkel excitons 

In spite of the above results, such a very visual way to calculate Shiva diagrams between N excitons and their 
decomposition in terms of exchanges between two excitons, encounter a major problem when dealing with the double 
exchange shown in Fig. 2(d): Indeed, Frenkel excitons Q[ and Q x are clearly made with the same electron- hole pair 
(see Fig. 2(e)); so that we should find ^Q^Q^Q^Qa as for Wannier excitons. In contrast, two Pauli scatterings between 
two Frenkel excitons give 



(of P ) A (p Qi) = ~pp E ^Qi+Qi.P'+P^P'+P.Qi+Cb 



p,p' p p 

i 



N S Q' x +Q' 2 ,Qi+Q2 ( 4 - 7 ) 

which is nothing but the Pauli scattering A ^q? q^J and thus definitely different from an identity. 

This problem must be associated to the fact that, while for Wannier excitons, the two possible ways to couple two 
electrons and two holes into two excitons leads to [12| 

4^] = -E A (- J i)^ ( 4 - 8 ) 

m n 

such a rearrangement is not possible for Frenkel excitons; so that a similar relation does not exist for them. Indeed, 
using Eq. (2.17), we do have 

44 = ^-E e^e^'aWM (4.9) 

n n' 

Due to the tight-binding approximation, we cannot make an exciton out of the electron-hole operator a^b^, except 
for n' — n, the product aj^a^fr^, then reducing to 0, according to Eq. (2.9). 

The difficulty to rewrite BqBq, in terms of two other Frenkel excitons, when compared to Wannier excitons, can 
be mathematically assigned to the fact that the free electron- hole pairs a k b\ \v) form a clean complete basis for one 
electron and one hole while this is not true for the trapped pairs a^b^ \v). As previously discussed, these states can be 
used as a basis in the tight-binding limit only, i.e., when the two lowest atomic (bound) states are the only relevant 
states of the problem. However, Eq. (4.8) as well as the identity between diagrams 2(d)-(f) for Wannier excitons, 
relies on an exact closure relation for all the intermediate electron-hole states on which the sums are taken while this 
closure relation using the two lowest atomic states is only approximate in the case of Frenkel excitons. 

This may lead us to think that, since there is one way only to construct two Frenkel excitons out of two electron-hole 
pairs, many-body effects coming from fermion exchanges between Frenkel excitons which result from their composite 
nature, must be somehow suppressed. The situation is actually far more subtle, as shown in the next sections. 

V. CALCULATION OF SOME RELEVANT SCALAR PRODUCTS 

As for Wannier excitons, physical effects dealing with Frenkel excitons ultimately read in terms of scalar products 
of Frenkel exciton states. Let us now calculate a few of them along the line we have used for Wannier excitons [l2l.[l9l|. 
Due to the structureless form of the Pauli scattering, these scalar products turn out to be far simpler than the ones for 
Wannier excitons. This makes the Frenkel exciton many-body physics induced by carrier exchanges also far simpler. 
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A. Scalar product of two-exciton states 

Using Eqs. (3.2,4), the scalar product of two-free-pair states reads as 

(v\ B„,B n ,Bi 2 Bl \v) = (v\ B n , {Bt 2 B n , - 6 n , „ 2 - D <n2 ] \v) 

ni^n f 2 712 ~^ ^n'-t ri2^n' 2 n± 2(5 n ^ n\^n' 2 n2^Tl\ n f 2 (5T) 

So that (v\ BnBniB^/B^ \v) = 1 — S n t n reduces to when n' — n as expected since a given site cannot be occupied by 
two electrons in the absence of degeneracy. This shows that while there are N s one-pair states B\ \v), the number of 
different two-pair states B^B^ \v) is N S (N S — l)/2 since to must be different from n, while B^B^ — B^B^. 
If we now turn to two-exciton states, the same procedure leads, using Eqs. (3.6, 9), to 

(v\ Bq^ Bq 2 B^ 2 -Bqj \v) = 5 Q[ Ql 5 Q > 2 q 2 + S Q[ Q2 5 Q , 2 Ql - Jj-Sq> 1+ q> 2 q 1+ q 2 (5.2) 

By noting that the last term of the above bracket is nothing but —A — A ^q, 1 , we see that this result 

which is shown in Fig. 3, is formally the same as the one for Wannier excitons, namely [l2| 

(v\ B m B n B\B) \v) = 6 mi S nj - A (^) + (to <-> n) (5.3) 

with (to, n, i, j) replaced by (Q' l5 Q' 2 , Qi, Q2). 

It follows from Eq. (5.2) that the norm of two Frenkel exciton states 

(v\ B^B^B^B^ \v)=l + S Ql Q2 - J- (5.4) 

differs from for Qi = Q2; so that there are different two-exciton states B^B^ \v) since Bq i Bq 2 — Bq 2 Bq 1 

while there only are N S (N S — l)/2 different two-pair states B^B^ \v). Consequently, as for Wannier excitons, the 
two- Frenkel exciton states form an overcomplete set [l6| . 

Note that, while the overcomplcteness of the two Wannier-exciton states readily follows from Eq. (5.3), we have 
shown that such a relation - which comes from the two ways to bind two free electrons and two free holes into two 
Wannier excitons - does not exist for Frenkel excitons due to the tight-binding approximation which is at the basis 
of the Frenkel exciton concept. Nevertheless, two-Frenkel-exciton states, like two- Wannier-exciton states, form an 
overcomplete set. 

B. Scalar product of three-exciton states 

If we now turn to the scalar product of three-exciton states shown in Fig. 4, we find using the commutators of Eqs. 
(3.6, 9) 

H b q> 3 Bq 2 B q/i B f Qi B^ B f Qs \v) = {S Q[ Ql 5 Q , 2 Q2 5 Q , 3 Qi + 5 perm.} 
2 



N f 



{^Qi+Q^ Qi+Q 2 ^Q^ Q 3 + 8 Perm.} 
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+ T^Qi+QJi+Qs Q1+Q2+Q3 ( 5 - 5 ) 

This result agrees with the standard rules [j| for calculating scalar products through Shiva diagrams shown in Fig. 
4, with the values of fermion exchanges between two and three Frenkel excitons, given in Eqs. (4.2, 4) and Eq. (4.5). 
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C. Scalar product of N identical exciton states 

As for Wannicr cxcitons, the normalization factor of N identical Frenkel excitons plays a crucial role in all physically 
relevant matrix elements involving N excitons, since most of them usually are in the same state. Let us write this 
normalization factor as [H, EH 

H B%B%> \v) = N\F N (5.6) 

while AH is the value this scalar product would have if Frenkel excitons were elementary bosons. We are going to 
show that, in contrast to Wannier excitons for which depends on the state (u, Q) in which the N excitons are, 
Fjy for Frenkel excitons depends on N but not on the exciton state considered, labelled by Q. This property is going 
to greatly simplify all results on many-body effects with Frenkel excitons. 
(i) Pedestrian way to calculate F N 

The pedestrian calculation of Fn relies on the commutator (3.10). While (v\BqBq \v) = 1 so that F\ — 1, it is 
easy to show that 

, t2 ,. -A 1 



(v\ B' Q B% \v) = 2 ( 1 - — 



,3 ot3. - ' - 1 W- 2 



(v\ B* Q B% \v) = 31 ^1 - W J [1 - W J (5.7) 
and so on... This leads us to conclude that Fjf should read as 

which makes Fn independent of the Frenkel exciton Q considered, a result somewhat surprising at first. 

(ii) More elaborate derivation 

A more elaborate way to calculate Fjy makes use of the recursion relation between the (v \ BqB^ \v) easy to obtain 
from Eq. (3. 13). This equation leads to 

(v\B$B}»\v) = (v\B»-iB Q Bl»\v) 
= {N-^j^)(v\B^B^\v) (5.9) 

So that 

F N =(l- F N -i (5.10) 

from which the expression of Fn given in Eq. (5.8) readily follows. 

(iii) Recursion relation between the Fjy's 

Even if Fn for Frenkel excitons has a nicely compact expression, it can be of interest to note that, due to Eq. 
(5.10), Fpf can be written as a sum of all the other F/v's according to 

JV-l (JV-l)(iV-2)_ 

Fn — Fff-i 77 — Fn^ 2 H 77^ Fn-3 — ■ ■ ■ (5-11) 

N s N* 

This series just is the analog of the one for Wannier excitons, namely (see Eq. (15) Ref. [HI]) 

F N = f2(-i r ^^^. anFN _ n (5.12) 

where the cr ra 's are the Shiva diagrams for fermion exchanges with n excitons Q on both sides. Indeed, as shown in 
Section IV, these Shiva diagrams in the case of Frenkel excitons reduce to l/iV" _1 (SQ' i+ ... + Q/ ii q 1 +...+q 71 , i.e., l/A^ 1 
when all the Q's are equal. 
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(iv) Density expansion 
From Eq. (5.10), we find 



F N -i N s ( ' f\ r 

-T^ = n s + i-n^ 1 + (n;) (5 - 13) 



for 1 <C N <C N s . This result is again fully consistent with Fm-i/Fm = 1 + 0(tj) as found for Wannier excitons fla ] 
where r\ — N (a x / L) with N being the number of excitons at hand, a x the exciton Bohr radius or relative motion 
wave function extension, L the sample size and D the space dimension. Indeed, the ratio of the sample volume L D 
divided by the Wannier exciton volume a® is nothing but the maximum number of excitons the sample can accept in 
the absence of spin degrees of freedom. This thus makes r\ for Wannier excitons equal to the ratio of the number of 
excitons at hand divided by the maximum possible number of excitons, this ratio being small in the dilute limit, i.e., 
when excitons exist. In the case of Frenkel excitons, it is clear that the maximum number of excitons has nothing to 
do with the Frenkel exciton relative motion extension because this extension reduces to zero in the tight binding limit. 
This maximum exciton number is just the number N s of atomic sites; so that the small dimensionless parameter for 
these excitons, associated to their density, must be identified with 

N 

(5.14) 

This expression leads us to reconsider the physical understanding of the parameter which controls carrier exchanges 
between composite bosons. While the spatial extension of the electron-hole relative motion wave function seemed at 
first to us a relevant length for the control of carrier exchanges, the expression of 77 we here find for Frenkel excitons 
shows that the physically relevant quantity for the low density expansion of carrier exchanges cannot be this length 
- otherwise 77 would reduce to zero in the tight-binding limit - but must be the maximum number of excitons we can 
put in the sample volume V . With this new understanding, it becomes easy to grasp why carrier exchanges are not 
totally suppressed in Frenkel exciton systems. 

If we now come back to the exact value of Fn given in Eq. (5.8), we can show that, as for Wannier excitons, it 
decreases exponentially with Nr). This is done by using the Stirling formula pi ~ p p e~ p y/2np in Eq. (5.8). For N s 
and N s — N both large, we find 

F N ~ e' N (l - V )- N ^~ 1 )(l - Ty)- 1 / 2 (5.15) 
Since a = e , we then get when all the Nrj's are small 

F N ~ e ~ Ar (^+T+-) ( 5 . 16 ) 
A similar exponential decrease is also found for Wannier excitons. 

VI. CLOSURE RELATION FOR FRENKEL EXCITONS 

We have shown that, due to the existence of a clean closure relation for the electrons and holes on which the 
Wannier excitons arc constructed, these composite boson excitons have a nicely compact closure relation which reads 

m 

1 = ^2 E 4 • ■ ■ B l \ v ) H Bit* • • • B h (6.1) 

the closure relation for elementary bosons being just the same with 1/ (JV!) replaced by 1/NL 

We are going to show that the situation for Frenkel excitons is rather different due to the tight-binding approximation 
underlying the Frenkel exciton concept. 

A. One electron-hole pair / one Frenkel exciton 

Starting from the closure relation for electrons in atomic site n 

I = Y, a n\v)(v\a n (6.2) 
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and the similar relation for holes, we can construct the closure relation for one-electron-hole pairs. It reads 

I = EE a ™ 5 ™'l w ) (v\K>a n 

n n' 

= E ^ \v) (v\ B„+EE a n b l \ v ) W & «' a « ( 6 - 3 ) 
Using Eqs. (2.18, 20), it is easy to show that the first sum can be rewritten as 

E B h> i«> (-1 b q w E e j(Q '- Q) - R " = E s q i«> (-i b q 

Q' Q S n Q 

while the second sum of Eq. (6.3) can be dropped out since it corresponds to projections over states way off by an 
energy e , their electrons and holes being in different sites. Consequently, the closure relation for one-Frenkel exciton 
states can be approximated by 

Q 

B. Two free pairs / two Frenkel excitons 

Let us now turn to two-pairs. The closure relation for two electrons reads 

1 = 2j E a «i a «2 \ v ) ( v \ a n 2 a ni (6.5) 
and similarly for holes. By rewriting this hole closure relation as 



/= ^E{ E blbl 2 \v)(v\ b< b n/i 



+ E b U b lj V ) Wn'M 
n^{ni,nj) 

and by splitting the sum over n' 2 in the same way, we find that the closure relation for two electron-hole pairs restricted 
to the degenerate subspace having an energy (e e + Sh — £o) can be written as 



J - \%l 2 E a «i fe «i a « 2 ^ 2 |«> (v\b n2 a n2 b ni a r , 



= \Y. B l B nM(v\B n2 B ni (6.6) 

If we now rewrite free pairs in terms of Frenkel excitons according to Eq. (2.18) and we use Eq. (2.20), the above 
closure relation leads to 

7 ^4 E B k B lM^\ B ^ B ^ ( 6 - 7 ) 

' Qi Q 2 

We see that this sum has the same prefactor 1/2! as the one for two elementary bosons. This equality between 
prefactors can be understood by noting that there is only one way to construct two Frenkel excitons out of two pairs 
belonging to the (s e + Sh — £o) subspace; so that these excitons behave as elementary bosons with respect to this 
projection operator. 
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C. Direct check 

It is possible to check the validity of the approximate closure relation for two free pairs, given in Eq. (6.6), by 
inserting it in front of \v). Indeed Eq. (5.1) for the scalar product of two-free-pair states readily gives 

\ E B nA l«> H ^4,44 \v) = B^Bl, \v) (6.8) 

for any m\ ^ m 2 , as expected for a closure relation, the state B}% \v) reducing to since the same site cannot 
accommodate two electrons. 

The check of the approximate closure relation for two Frenkel excitons, given in Eq. (6.7), is less trivial. Indeed, 
if we insert Eq. (6.7) in front of £>Pi-^p 2 \v), we find, using Eq. (5.2) for the scalar product of two- Frenkel exciton 
states, 



\ E B k B k H H^q^q^pA \v) = 4 A l«> 



2 
QiQ 



Qi Q 2 



In order to show that the above sum reduces to zero as necessary to prove that Eq. (6.7) is indeed a closure relation, 
we use Eq. (3.11) for A (g* We find 

E pO« = 1v-E4 1+ q4 2 -qH (6-10) 



Qi Q 2 N " Q 

We then rewrite exciton operators in terms of free pairs according to Eq. (2.17) and we use Eq. (2.19). This leads to 

E4 1+q 4 2 -q = E ^ p ^» 1 ^-^)Bt iB t 2 ^ e ^.(H„ 1 -H„ 2 ) 

Q nin 2 Q 

= N s Y / e i{Fl - P2) - R ^Bt 2 i (6.11) 
So that since, B 2 = 0, due to Eq. (2.9), we end with 

E x {t p0« = ° ( 6 - 12 ) 

Qi Q 2 

This result must be physically related to the fact that there is only one way to write two Frenkel excitons out of two 
pairs: If we try to exchange the carriers of two Frenkel excitons as in the sum (6.12), we get zero. 

D. Closure relation for N Frenkel excitons 

A similar procedure performed for N pairs leads to an approximate closure relation for N Frenkel excitons which 
reads 

7 - E B Qi ■ ■ ■ B Qn \ v ) H B Q« • • • B Q> ( 6 - 13 ) 
with again the same I/AH prefactor as the one for N elementary bosons. 
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VII. CREATION POTENTIAL FOR FRENKEL EXCITONS 

Frenkel excitons feel each other through exchanges as described in the previous sections. This "interaction", 
generated by the Pauli exclusion principle, is of course unusual in the sense that it gives rise to dimensionless Pauli 
scatterings as shown in Section III. Frenkel excitons also interact through the Coulomb potential which exists between 
their fermionic components. As for Wannier excitons [l3l |. there is no clean way to write Coulomb interactions between 
two excitons as a potential due to the composite nature of the particles. Here also, we are going to overcome this 
difficulty by introducing a set of " creation potentials" . From them, it becomes formally possible to derive the 
interaction scatterings of two Frenkel excitons, as we shown in the next section. 



Definition 



The creation potential Vq of the Frenkel exciton Q is defined through 



— EaB, 



V, 



Q 



(7.1) 



where Hx = H 



(0) 
X 



Vcoui is the interacting exciton Hamiltonian defined in Eq. (2.24) while 



H, 



eh 



Sx with 



H e h, Sx and V cou i given by Eqs. (2.3, 22, 25). To get we have to calculate the commutator of the Frenkel exciton 
creation operator Bq with each of the three terms of Hx- However, since (iJ^ 0) - £ Q ) B+j \v) = 0, the EqB^ part 



in the RHS of Eq. (7.1) must come from 



o-(o) R t 



. This leads us to split Eq. (7.1) as 



tl x ,±S, 



Q 



— EaB, 



(7.2) 



Vcoui i Bq 



(7.3) 



The Frenkel exciton creation potential is thus going to appear as 



(7.4) 



these two parts being conceptually quite different: Wq associated to V cou i comes from direct Coulomb processes 

between sites. It has similarity with the Wannier exciton creation potential. In contrast, Sq, which originates from 
Sx, essentially comes from indirect Coulomb processes; these processes, crucial for Frenkel excitons as they are 
responsible for the excitation transfer, are usually neglected for Wannier excitons. They precisely correspond to the 
so-called "electron- hole exchange". 



B. Free exciton creation potential 



The free exciton Hamiltonian contains an electron - hole kinetic energy term H e h defined in Eq. (2.3). Its 
commutator with £>q readily gives 



H e h, Bq 



£ e ^2 a\ L a n +e h ^2 b li h n, ^2 



7K 



'■II i i > i 



= (e e + E h )B^ (7.5) 

H^ also contains Coulomb contributions through Sx ■ Although Sx reads as a one-body diagonal potential in the 
Frenkel exciton subspace, it can produce interactions between two excitons through the carrier exchanges the exciton 
composite nature allows - an idea not trivial at first. 
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By using Eq. (3.6), it is straightforward to see that 
for Frenkcl exciton as 



Sx,B, 



Q 



reads in terms of the deviation-from-boson operator 



Sx,B Q 



Q' 



ECQ'< (^Q'Q-^Q'Q) 



Q' 



= CQ^-ECQ' 5 Q^Q'Q (7-6) 
Q' 

When combined with Eq. (7.5), we readily get Eq. (7.2), with Sq given by 

S Q = - E (Q' B h' D Q'Q - - E B Q'tQ' ( A Q-Q' + A Q-Q') ( 7 ' 7 ) 
Q' Q' 

Cq appearing in Sq is defined in Eq. (2.23). It has a constant contribution (— e ) which comes from the neutrality 
potential Vi n t ra induced by intrasite Coulomb interactions. It also has a contribution Vq which comes from electron- 
hole exchange interactions between sites which allow to delocalize the excitation. As shown below, the interaction 
scatterings generated by these two contributions have rather different properties; this is why it is physically relevant 
to split Sq as 

Sq=^q+ t q (7-8) 
where iVX and 7q read as Sq in Eq. (7.7), with (q> replaced by (— e ) and Vq' respectively. 



C. Interacting exciton creation potential 



We now turn to the creation potential induced by direct processes through V cou i. This interaction contains three 
parts. The electron-electron contribution gives 



V ee , B Q 



iQ.R„ 



^ E ^»i-R. J (il) fl Iifll'^flni. E 77fi= Q " 6 " 



(7.9) 



By using [a^a^a^a^, a^b^] = a\ n b\ n {8 nim a\ l2 a n2 + 6 n2Tn al i a ni ) and by writing a^b^ in terms of excitons 
through Eq. (3.10), we find, since V R ({ \) = V- R (\ {) 



^e, B Q ] =E S Q'E aWe^-«^ E ^^^^ G I) (7-10) 

Q' n m=jtn 



The last sum does not depend on n due to the translational invariance of the system; so that the above equation can 
be rewritten as 



V ee , B Q ]=E«-Q' A 
Q' 



(e) 

Q Q' 



(7.11) 



where Wq 6 ' is very similar to the Q component of the exciton energy Vq given in Eq. (2.15), except that it reads in 
terms of direct Coulomb interactions between two excited atomic levels v = 1. Its precise value is given by 



< e) = E e ~ iQR ^ i) 

In the same way, the hole-hole contribution of V cou i leads to 

Vhh, B Q ]=Y, B h^ Q h %A^_ Q , 
Q' 



(7.12) 



(7.13) 
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where W Q hh) reads as W Q ee) with V R (\ \) replaced by V R (g g). 

If we now turn to electron-hole direct interaction, its commutator with Bq gives 



Tr(dir) R f 
V eh ' n Q 



E^R--R»'(S5)W'«n'6n. E 



u m u m 



(7.14) 



To calculate it, we use 



bla^an'bn, a^b^ = aJ„&J rl ((5„ m (5„'„ l + 5 nm a], ,a„/ + (Jn'm&Jj&n)- The first term of this commu- 
tator gives for n ^ n'. In the two other terms, we rewrite aj^fr^ in terms of excitons according to Eq. (2.18) and 
we make appearing e 4 ^ _Q )-( R ™~ R ™) in the remaining sum to possibly use the translational invariance of the system. 
This ultimately leads to 



v (dir) R f 



E B l E KVS^Q* + <- } Q' A ^Q') 

Q' n 



(7.15) 



where J and read as W^ eeJ given in Eq. (7.12) with Vr ({ \) replaced by Vr (J J) and V R (? ?) respectively. 

By collecting these three contributions, we find that Wq 



VcquiBq 



reads as 



K =E S Q' ( W Q-Q' A Q-Q' + W Q-Q' A Q-Q') 



Q' 

where the prefactors are such that W Q e) = wj e) - W§ h) while W Q h) = W Q hh) - W Q he) ; so that 



J) = ^ e -iQ.R [yk ( 1 1) _ Fr ( l 1 }] 

R/0 



(7.16) 



(7.17) 



^ - E *" <Q R 

R/0 



(7.18) 



D. Discussion 



It is of importance to note that, while the structure of the two parts Sq and Wq of the Frenkel exciton creation 
potential Vq are rather similar, the prefactors of Bq, in Wq which come from direct Coulomb processes between 
electrons and holes of different sites, depend on the momentum transfer Q — Q' only, while this is not true for the 
prefactor of Bq, in Sq which comes from " electron- hole exchange", i.e., indirect Coulomb processes between atomic 
levels. This result, which may appear as very strange at first, is linked to the fact that these indirect Coulomb processes 
do not lead to a momentum transfer: Indeed, the potential Sx they produce is a one-body diagonal operator in the 
exciton subspace, as seen from Eq. (2.22). 

Let us recall that the Coulomb exchange processes producing Sq in which one electron is excited while the other 
returns to its ground state, are usually neglected in the case of Wannier excitons. Being much smaller than the direct 
scatterings due to the orthogonality of the ground and excited state wave functions, they just produce the small 
energy splitting which exists between dark and bright excitons, when the spin degrees of freedom of these Wannier 
excitons are included. In the case of Frenkel excitons, the prefactors in Sq are also small compared to the ones coming 
from direct processes since through Vr (5 J), they also contain overlaps of different atomic states v = and v = 1 - 

which are orthogonal - while the prefactors of Wq contain Vr Jf), i.e., overlaps of identical atomic states, v or v' . 
However, these indirect Coulomb processes are crucial for Frenkel excitons since, in the tight binding limit, they are 
the only processes allowing to delocalize the excitation, as necessary to produce the exciton. 
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VIII. INTERACTION SCATTERINGS 



Definition 



The interaction scatterings of the two Frenkel excitons (Qi, Q2) are defined through the creation potential Vq as 



E 

Qi Q' 2 



€ (qI Qi) B Qi B Q 2 



(8.1) 



Before going further, let us recall that this equation, when used for Wannier excitons, does not unambiguously 
define [20] the interaction scattering £ between (i, j) and (m, n) states. Indeed, due to the two possible ways 

to write two excitons leading to Eq. (4.8), it is always possible to replace £ fm^i) m (8-1) by mmus the " m " 
exchange interaction scattering defined as 



C n ( n j .) = Va ( n q )e ( q 3 

p q 



.2) 



or even by |a£ (™ ^ - 6f " (™ ^ J J with a + b = 1. 

In contrast, due to the fact that there is only one way to construct two Frenkel excitons out of two electron-hole 
pairs, an equation like Eq. (4.8) does not exist for Frenkel excitons. This leads us to think that the value of the 

interaction scattering £ ^q? q^ obtained through Eq. (8.1) must be unique. As a direct consequence, the "in" 

exchange interaction scattering for Frenkel excitons must be equal to zero. A way to show it, is to insert the closure 
relation (6.7) for two Frenkel exciton states in front of RHS of Eq. (8.1) acting on vacuum. We find 



Qi Q 2 



Pi p 



where the scattering appearing in the sum is given by 

?(p^H E ^(q?q:)h^^«i-) 



(8.3) 



(8.4) 



Qi Q 2 



If we now use Eq. (5.2) for the scalar product of two Frenkel exciton states and we symmetrize the interaction 
scattering, £ (^f Qi) = £ (q' 1 Qi) > as arwavs possible since £?q, B^, = B^, B^, , we find that £ q^ is just 



7(P 2 Q 2 \ _ c (P 2 Q 2 \ _ Hn (P 2 Q 2 



■5) 



where £ m q^ is the "in" exchange interaction scattering of two Frenkel excitons, defined in the same way as for 
Wannier excitons, namely 



tin fPz Qz\ _ \ (Pi Q 2 \ f [Q 2 Q2 

? {p, Q X J - A \Pi Q'J ? {Q[ Q 1 

Qi Q 2 



■6) 



By inserting £ given in Eq. (8.5) into Eq. (8.3), we are led to conclude that £ m = 0, or if £ m has a non zero 

part, its contribution to the sum (8.3) must give zero. 

Let us now see how this conclusion, based on very general properties of Frenkel excitons, can be recovered through 
hard algebra. 



B. Calculation of the interaction scattering 



To calculate the interaction scattering £ ^q? q^ defined in Eq. (8.1), with the creation potential Vq given in Eq. 
(7.4), we first note that 



A (e) £ f 



E 



Ve lQ R "V fet 
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(8.7) 



So that, by writing free electron-hole operators in terms of excitons according to Eq. (2.18), and by using Eq. (2.20), 
we end with 



This readily shows that 



A (e) R f 



A(h) „t 



Q2 



~~ M X!^i jB Qi B Qi+Q2-Q' 1 



Qi 



(8. 



(8.9) 



w L B k 



N, 



Qi 



(8.10) 



where Cq is defined in Eq. (2.23) while Wq is equal to Wq + W%'; so that 

wq = E e " iQ R G \) + (g °o) v R (I I) V R (? ?)] 

R#0 



(8.11) 



It contains all possible direct interactions between electrons and holes in different sites. This quantity is real, since 
Vr, ") = V- R (" while (J J) = F-r (5 J). For the same reason, it is an even function of Q. 

From Eqs. (7.4) and (7.8), we are lead to split the interaction scattering defined in Eq. (8.1) into three conceptually 
different terms 



c (Q' 2 QA _ t (Q' 2 QA _ t (Q 2 Q 2 \ , (Q' 2 Q 2 
? VQi Qi / q,coul \Qi Qi ) ?trans \Qi Qi J q,neut \Qi Qi 



(8.12) 



The first contribution which reads 



t (Q' 2 Q2 

C,coul lq; q 1 



- ^T^Qi-Qi^Qi+Qi Q1+Q2 



(8.13) 



only depends on the momentum transfer Q x — Q' x produced by the scattering. Through Wq, this contribution comes 
from all direct Coulomb processes between the excited atomic level v = 1 and the ground state level v = of the 
different sites. This part of the interaction scattering which is represented by the diagram of Fig. 5(a), is very similar 
to the direct Coulomb scattering of two Wannier excitons. 

The two other contributions to the interaction scattering are conceptually new. By symmetrizing Eq. (8.9), it is 
possible to write them as 



£,transf (q/ q") - ^f^Qi + V Q 2 ^Qi +Q2 >Qi +Q2 



5.14) 



t (Q' 2 Q 2 \ _ 

Z,neut ^q/ Ql J - 



2e 0x 

-J^-OQ' 1 +Q' 2 ,Qi+Q 2 



(8.15) 



These scatterings may appear as very strange at first since they do not depend on the momentum transfer but on the 
momenta of the two "out" excitons Q[ and Q' 2 for £transf, while £ neu t is completely independent of the momenta of 
the Frenkel excitons involved in the scattering - which is even stranger. 

Actually £„ ettt does not play a role in the scattering of two Frenkel excitons. Indeed, when inserted in Eq. (8.1), 
its contribution gives zero 



Ep Q2A R t R t _ n 

t,neut y Qli Ql j n^n^ - u 



(8.16) 
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This can be readily shown by noting that 

Ueut (q? Qi)= 2£ A (q; q') ( 8 - 17 ) 

and by using Eq. (6.12). Consequently £ in Eq. (8.1) can as well be replaced by 

? {Q[ Qi) ~ ?cou/ \Q[ Qi ) ?trans \Q[ Qi J \o.io) 
If we now turn to ( tran5 j, given in Eq. (8.14), we can note that, due to Eq. (3.11), it also reads 

Ztransf £) = (Vq/ + Vq/)A (g J) (8.19) 

While its contribution in the sum (8.1) definitely differs from zero, this expression of £,transf leads us to see it as 
a "transfer assisted exchange" (see Fig. 5b). And indeed, it appears in the interaction scattering £ or £ with the 
minus sign as for any exchange process. To better grasp this term, we can note that, although the potential Sx from 
which it originates, is a one-body operator in the Frenkel exciton subspace, it can lead to a scattering between two 
excitons thanks to carrier exchanges allowed by the exciton composite nature. Such a scattering Straus /, which does 
not exist for Wannier excitons when electron- hole exchange, i.e., valence-conduction Coulomb process, is neglected, 
is quite specific to Frenkel excitons: Indeed, besides direct Coulomb processes which exist for both excitons, Frenkel 
excitons in addition have indirect Coulomb processes which provide other energy-like quantities from which this novel 
energy-like scattering can be constructed, by mixing them with carrier exchanges. 

It can be of interest to relate this ^transf scattering to the novel scattering we have identified between two polaritons 
[HI - we called " photon assisted exchange" . For polaritons, the Pauli scattering induced by carrier exchanges between 
the exciton components, is mixed with the Rabi coupling between photon and exciton, which also is a relevant energy- 
like quantity appearing in the one-body part of the polariton hamiltonian. As for this photon assisted exchange, the 
strange asymmetry between "in" and "out" states of assisted exchange scatterings is actually cured in real physical 
processes in which energy is conserved since the scatterings with interaction before or after the exchange are then 
equal, as shown in Section IX. 



C. Exchange interaction scatterings 

The exchange Coulomb scatterings which are succession of carrier exchange before or after carrier interaction, play 
a key role for Wannier excitons. They are the dominant scatterings for excitons with close-to-zero momenta. We are 
going to show that in the case of Frenkel excitons these exchange interaction scatterings reduce to zero 

rg £) = E*($ $) = (8 - 20) 
PiP 2 

By using Eq. (3.11) for the Pauli scattering, we readily find that the "in" exchange scattering associated to direct 
Coulomb processes given in Eq. (8.13) reduces to zero 

Cuz (2? 20 = j^Qi+Qi.Qi+Q* E W Q = ( 8 - 21 ) 

Q 

due to Eq. (2.19), namely ^2q e *Q- R — N s 6-r.o since the R = term are also missing for the sums which enter Wq 
(see Eq. (8.11)). 

In the same way, the "in" exchange scattering associated to the transfer part of the interaction given in Eq. (8.14) 
reduces to zero 

Ctransf [q{ 2l) = Jj2 ^Qi+Qs.Ql+Q^ EO^Ql+Q + V Q 2 -q) = ( 8 - 22 ) 

Q 

since the R = term is also missing from the sum entering Vq (see Eq. (2.15)). 

We thus conclude that, as obtained from very general arguments at the beginning of this section, the "in" exchange 

scattering for the first two parts of £ (J^ q\)i namely (^f ^\ and £J"ans/ (2' 2i) rea ^ uce to zer0 while the 
third part ^ neu t gives zero when inserted in the relevant sum (see Eq. (8.16)). This physically means that the energy- 
like part of the Frenkel exciton scattering is only controlled by direct Coulomb processes and by transfer assisted 
exchange, a result not obvious at first. 
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IX. FRENKEL EXCITON HAMILTONIAN IN THE TWO-EXCITON SUBSPACE 

A. H x acting on two excitons 

In order to better grasp the interplay between the various scatterings appearing in the preceding sections, let us 
calculate the matrix element of the interacting exciton Hamiltonian in the two-exciton subspace. Using Eq. (7.1), we 
find that Hx acting on two Frcnkel excitons gives 

HxB^B^ \v) = (B^Hx+E^B^+V^) B^ \v) (9.1) 

Due to Eq. (8.1), this also reads 

(Hx-E Ql -E Qa )B^ Qa \v) 

= « l«> = E € (£ Qi) « l«> (9-2) 
Pi p 2 

where £ can be possibly replaced by £ defined in Eq. (8.18), due to Eq. (8.16). 

B. Matrix element of Hx in the two-exciton subspace 

We now turn to the matrix element of the interacting exciton Hamiltonian Hx in the two exciton subspace. With 
Hx acting on the right, wc find 

(v\ B^B^HxB^B^ \v) = (E Ql + E Q2 ) (v\ B^B^B^B^ \v) 

+ E H s q^q;« \v) £ (£ t) (9-3) 

Pi p 2 

If we now use the scalar product of two-exciton states given in Eq. (5.2), the first term of Eq. (9.3) reduces to 

(E Ql + Eq 2 ) ^<5q' 1 q 1 <5q^q 2 + ^QiCh^QjQi - ^ <5 Qi+Q 2! Qi+Q2^) ( 9 - 4 ) 

To calculate the second term, we use the same scalar product and take into account that the exchange scatterings 
reduce to zero according to Eq. (8.20). This second term then reads 

^{q! qI)+KqI 20 = ]^( W Qi-^ +W ^-^- 2V Qi- 2V ^« y Q' 1 +Qi.Qi+Q a (9-5) 

C. On the asymmetry of transfer assisted exchange scatterings 

In order to better accept the asymmetry between " in" and " out" excitons appearing in the transfer assisted exchange 
scattering, let us now show that these "in" and "out" scatterings are equal when energy is conserved. For that, we 
use the above equations to calculate the same matrix element (v\ Bq^Bq^HxB^B^ \v) but through 

{v\B Q2 B Ql HxB^,B^, 2 \v)* 
+ ^ (Wqx-QJ + Wq 2 -Q 2 - 2V Q 2 - 2 ^Qi) <Wq 2 ,Q'i+Q 2 



23 



By noting that Eq is real as well as Vq (see Eq. (2.15)), while Wq = Wq = W-q (see Eq. (8.11)), and by comparing 
the two expressions of (v\ Bq^Bq^HxB^B^ \v), we readily see that 

(Vq; + V Q > 2 - V Ql - Vq 2 ) <5 Q ' 1+ q^q 1+ q 2 

= (Eq^ + E Q > 2 - Eq 1 - Eq 2 ) 5 Q > i+ Q> 2tQl+ Q 2 (9.7) 



which also reads 



t (Q' 2 Ch\ _t (Q2 Q' 2 \ 

Straus f ^q/ q 1 J C,transf q^ J 



= {E Q[ + E Q , 2 Eq 1 Eq 2 ) A g? J) (9.8) 

This proves that the asymmetry between " in" and " out" states in the transfer term is removed for energy conserving 
processes which are the ones controlling the physically relevant processes between two Frenkel excitons. A similar 
result is found for the photon assisted exchange scattering of two polaritons. 

D. Hamiltonian expectation value 

From Eqs. (5.2) and (9.3-5), we get the Hamiltonian expectation value for two Frenkel excitons as 

(v\B Ql B Q2 H x Bl Qi Bl Q3 \v) 



(vlB^B^B^B^Jv) 



= Eq, + E Q2 



1 W + >V Q2 _ Ql -2V Ql -2V Q2 
+ N S 1 + SqiQ2 -± 

So that the Hamiltonian expectation value for two identical excitons Qi = Q2 = Q reduces to 

(H X ) 2 = 2E Q + (9.10) 

Its interaction part contains a W contribution which comes from direct Coulomb processes between sites. It also 
contains an indirect "transfer" contribution induced by the possible fermion exchanges between composite boson 
excitons. This term which comes from exchange between two excitons mixed with excitation transfer has a nature 
quite different from the one of the direct Coulomb term Wo- As Eq, it depends on the momentum Q of the exciton 
considered. We sec that the interaction part of the Hamiltonian expectation value of two excitons reduces to zero in 
the large sample limit, i.e., for N s — ► 00. This is physically reasonable since the interaction of just two excitons in 
a huge volume cannot change their energy very much. However, as in the case of Wannier excitons, this interaction 
part is going to appear with a N(N — l)/2 prefactor when considering N excitons instead of two, due to the number 
of ways to choose 2 excitons among N. The precise calculation of the Hamiltonian expectation value for N Frenkel 
excitons is out of the scope of the present paper. 

E. Comparison with elementary bosons 

Since many people are still tempted to replace excitons by elementary bosons with a set of " appropriate" effective 
scatterings dressed by exchanges, let us again show that, as for Wannier excitons, it is not possible to construct such 
an effective bosonic Hamiltonian for Frenkel excitons, due to possible carrier exchanges between these excitons. 

If an effective Hamiltonian were to exist, it would read 

Hx=Y, E F bIb p + \ Yy (p? £) Bp'X, B P2 B Pl (9.11) 
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with 



<5p<P for boson excitons; so that H x acting on two-exciton states would give 



*MqX 2 l«) = (^Qi +^Q 3 )^ t Q 1 ^Q 2 \v) 



V 



fQ'2 QA\ pt -nt 



(9.12) 



To get the effective scatterings V f q? q* 



we can think of comparing this result with Hx acting on two Frenkel 

excitons. By symmetrizing Eq. (9.2) for Bq, Bq, = Bq, £?q, , this leads us to take V ^q? q^ as £ (^J or better 

as £ ^q? q^J in order to avoid the spurious constant term (— 2eq) in the effective scattering. However, this result 
is fully unacceptable since the (Vq' + Vq^ ) part of Straus / induces a non-hermiticity in the effective Hamiltonian. 
Indeed, for Hx to be hermitian, we must have 



V 



p' 2 p, 
P', P: 



V 



(9.13) 



while Straus f does not have this property. Consequently, as for Wannier excitons [20(, it is not possible to construct 
an effective bosonic Hamiltonian for Frenkel excitons with the structure given by Eq. (9.11), due to the presence 
of exchange processes in the matrix element. However, we would like to notice that in Refs. [2l|, [2^, [H| the 
fermionic Hamiltonian was represented exactly as an infinite series expansion in creation and destruction operators 
for excitations on atomic sites and these operators obey bosonic commutation rules. Note that, within such a scheme, 
it is rather complicated technically to treat processes involving large number of excitons, as explained in a more detail 
in Section XI. 



X. DISCUSSION 



Since Frenkel excitons are made of highly localized electron-hole pairs, we can be led to think that carrier exchanges 
between Frenkel excitons should be strongly reduced; so that their exciton composite nature should have very little 
consequences on their many-body physics. To check this basic idea, we are going to make a precise comparison 
between Frenkel excitons and elementary bosons. 

The fact that these excitons are made with electrons and holes localized on the same site, should also make them 
rather different from Wannier excitons constructed on delocalized conduction electrons and valence holes. A precise 
comparison of the basic properties of Frenkel and Wannier excitons is thus of interest to better grasp the relevant 
parameters which control the many-body physics of these composite bosons. 



A. Comparison with elementary bosons 

Let us first compare Frenkel excitons with elementary bosons. 

Frenkel excitons are made of indistinguishable electrons and holes. They thus are composite particles. However, 
due to the fact that the wave functions of these carriers are the ones of highly localized atomic states on the same 
ion site, there is one way only to construct two excitons out of two electron-hole pairs. This leads to make Frenkel 
excitons appear as elementary particles. To support this idea, we can mention their closure relation which has the 
same (1/AT!) prefactor as the one for elementary quantum particles. 

In spite of these arguments, Frenkel excitons differ from elementary bosons by various means. A very important 

aspect of their composite nature is the existence of Pauli scatterings A ^q? q^J between two Frenkel excitons induced 

by carrier exchanges. These dimensionless scatterings in particular appear in the scalar product of two Frenkel exciton 
states, which does not reduce to a set of Kronecker symbols, as for elementary bosons. The possible carrier exchanges 
between excitons also produce a normalization factor for N identical excitons which is not Nl as for elementary bosons 
but NIFff where Fx is not a small correction but decreases exponentially with N. Finally, when mixed with the 
indirect Coulomb processes responsible for the Frenkel exciton excitation transfer between sites, these dimensionless 
Pauli scatterings give rise to an energy-like "transfer assisted exchange scatterings" which cannot be properly included 
in a bosonic effective Hamiltonian for Frenkel excitons without dropping its required hermiticity. 

These Frenkel excitons thus are somewhat awkward because they partly behave as elementary particles and partly 
as composite particles. This makes their handling through intuitive arguments rather dangerous. This is why the 
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hard algebra based on four commutators proposed in the present many-body theory, appears as highly valuable to 
overcome any possible uncertainty. 

B. Comparison with Wannier excitons 

Let us now compare Wannier and Frcnkel excitons. 

Wannier excitons are constructed on two fully delocalized particles, the conduction electron and the valence hole. 
The problem is then to make out of these delocalized fermions a bound exciton. This is done through intraband 
Coulomb processes. In contrast, Frenkel excitons are constructed on two fully localized particles, the first atomic 
excited state of a given site and the absence of atomic ground state electron on the same site. The problem then is to 
delocalize this highly bound excitation to produce a delocalized excitation. This is done through interband Coulomb 
processes. 

Wannier excitons are characterized by a center-of-mass momentum and a relative motion index, the electrons and 
holes on which these excitons are constructed being two fully delocalized particles. In contrast, Frcnkel excitons arc 
only characterized by a center-of-mass momentum since their electrons and holes are not free but forced to be on the 
same site. Due to this unique degree of freedom, N Frenkel excitons have the I/AH prefactor in their closure relation 
characteristic for N elementary particles, while the closure relation for Wannier excitons has a prefactor (1/AH) 2 which 
is the signature of the two degrees of freedom of these excitons. 

Being both composite particles, they can interact by fermion exchanges through 2x2 Pauli scatterings. These 
fcrmion exchanges make the scalar product of N identical excitons differ from their elementary boson value N\ by a 
factor Fn, which, for both excitons, decreases exponentially with Nij. The small dimensionless parameter r\ associated 
to density N/L D which controls these fermion exchanges, appears as N (a x /L) D for Wannier excitons, while it appears 
as N/N s for Frenkel excitons. This proves that the key parameter for the many-body physics of excitons is not linked 
to the extension of the electron-hole relative motion wave function as first thought - otherwise r\ would reduce to 
zero in the case of Frenkel excitons. In contrast, it must be physically understood as the number excitons at hand 
divided by the total number of excitons the sample can accommodate, this number being the number of ion sites N s 
in the case of Frenkel excitons while it is the ratio (L/a x ) D of the sample volume to the exciton volume in the case 
of Wannier excitons. 

The Pauli scatterings for Wannier and Frenkel excitons are represented by the same Shiva diagram for carrier 
exchange between two excitons. This diagram, as well as the ones for carrier exchange between more than two 
excitons, can be decomposed in terms of 2 x 2 Pauli scatterings and calculated along the same intuitive line as the one 
we have established for Wannier excitons. The only "little" problem arises with two consecutive exchanges: Instead 
of getting an identity, as visually seen, two consecutive exchanges between two Frenkel excitons is found to reduce to 
an exchange. As a direct consequence, it is not possible to rewrite two Frcnkel excitons as a sum of two other Frenkel 
excitons, as in the case of Wannier excitons. This awkward result on two consecutive exchanges puts some shade on 
the very visual handling of Shiva diagrams in the case of Frenkel excitons. This is why it seems wise to support all 
visual results on Shiva diagrams by a hard algebra of the Frenkel exciton matrix element they represent, using the 
four key commutators of the present many-body theory. 

The possible difficulties with Frcnkel exciton Shiva diagrams, when compared with Wannier excitons, can be traced 
back to the fact that the states on which Frenkel excitons are constructed, do not form a clean complete basis for 
all electron-hole pair states due to the tight binding approximation on which these excitons are based. This is why, 
when summations over intermediate states have to be performed, they not always give the correct results. To predict 
when these results are incorrect is however not so easy. Thanks to the full proof procedure we have constructed, it is 
however possible to reach the correct answer through a rather simple algebra based on commutator manipulations. 

Besides fermion exchanges, Frenkel and Wannier excitons also fell each other through Coulomb interaction. In 
addition to direct processes between carriers, indirect processes are also kept for Frenkel excitons since, in the tight- 
binding limit for electron and hole states, these are the only processes allowing to delocalize the excitation, i.e., to 
produce the exciton. Although appearing as a one-body operator in the exciton subspace, these indirect "transfer" 
processes give rise to a novel 2x2 energy- like scattering, when mixed with carrier exchange. These "transfer assisted 
exchange" scatterings which do not exist for Wannier excitons, are specific to the Frenkel exciton many-body physics. 
They in particular appear in the matrix element of the Hamiltonian in the two exciton subspace as well as in the 
Hamiltonian expectation value. 
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XI. STATE OF THE ART 

Some many-body properties of Frenkel excitons have already been studied in the past, see Refs. [2ll I22I I23I l24l. I25I 

M- _ __ 

(i) Agranovich and coworkers [2l], [22, [23| , as well as Davydov in Ref . [24] , start with a system Hamiltonian written 
in terms of creation operators a\ n and a\ n for ground and first excited states of isolated molecules on site n which is 
quite similar to our Hamiltonian with ao n essentially replaced by w n . Their fermionic Hamiltonian is then rewritten in 
terms of excitation operators B^, i.e., products like a\ n ao n - To do so, it is invoked [2l[ that the sum of the occupation 
numbers for the ground and first excited states of each molecule is equal to 1. This condition between scalars is then 
replaced by the same condition between operators, as seen in Eq. (4.9) of Ref. [2l|, namely 

7Vo„ = 1 - N ln (11.1) 

where Nf n — oX n cifn is the occupation number operator of the state / = (0, 1) in the site n. Eq. (11.1) is valid only 
when applied to states in which the sum of the occupation numbers for the ground and first excited states of each 
molecule is equal to 1. This equation (11.1), when implemented to the theory significantly facilitates a final compact 
form for the Hamiltonian. As seen from Refs. [U [22j , this procedure has however been used for two-exciton problems 
only. The generalization [23| of this approach to an arbitrary number of excitons turns out to be rather complicated 
technically, in contrast to our method which is quite straightforward for any number of Frenkel excitons, thanks to 
Eqs. (3.12, 13). 

(ii) Mukamel and coworkers [ID, [2(| H3, HI] have followed a completely different approach based on the idea that, for 
problems dealing with a fixed number N of electron-hole pairs, it is always possible to write the system Hamiltonian H 
as well as the commutators between electron- hole operators [B m , B^\ exactly, in terms of infinite series of electron-hole 
operators B^... B^ p B m B mp , namely 

00 

p=l 

H (v) = h^(n 1 ...n p ;m 1 ...m p )Bl l ...Bi p B mi ...B mp (11.2) 

{n}{rn} 

and similarly for \B m , . In order to describe a system of N interacting composite excitons, we then have to keep 
terms in these series for the Hamiltonian and the pair commutator, which contain products up to the p = N and 
p = N — 1, respectively. The constant prefactors in front of all these operators are obtained through the projection 
of the initial fermionic Hamiltonian and the commutator [B m ,i?^] on the A-pair subspace. This method which is 
formally exact, turns out to be technically quite heavy since the series become more and more complicated with the 
increase of p. Actually, the authors of Refs. [25|, [26|, [2?], HH have derived these series for two-exciton problems only, 
and they have used them to calculate successfully the third order susceptibility x' 3 ' ■ 

Very recently, we have reconsidered the idea of using infinite series of operators to describe interacting composite 
excitons (29j . Instead of using free pair operators, like Mukamel and coworkers, we have used correlated pair operators, 
i.e., excitons which are the relevant operators for such problems. By making use of the exciton closure relation, we 
have previously derived, we have been able to write the two infinite series explicitly, through rather light calculations. 
These series read in terms of Pauli and interaction scatterings. However, the proper handling of these two infinite 
series for A-body problems turn out to be far more cumbersome than the many-body theories for Wannier or Frenkel 
excitons we have constructed, based on just four commutators. 

XII. CONCLUSION 

In this paper, we propose a new many-body theory for Frenkel excitons in which the composite nature of these 
particles is treated exactly, excitons being never bosonized. Starting from the expression of the Frenkel exciton 
creation operator Bq written in terms of electron and hole creation operators, we derive the commutation rules for 
these Frenkel excitons. They clearly differ from the ones of elementary bosons due to carrier exchanges induced by 
the exciton composite nature. These carrier exchanges give rise to "Pauli scatterings" between excitons, similar to the 
one we found for Wannier excitons. They are represented by the same Shiva diagrams calculated along the same line. 
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Due to the fact that Frenkel excitons are only characterized by a center-of-mass momentum Q, their Pauli scatterings 
take an extremely simple form 



A 



(2i Qi) = jv^Qi+Qi. Qi+Q^ I 12 - 1 ) 



where N s is the number of ion sites, the many-body physics of N Frenkel excitons being controlled by the small 
dimensionless parameter 

V = N/N s (12.2) 

This parameter is the analog of parameter r\ = N (a x /L) D for Wannier excitons since N s like (L/a x ) D is the maximum 
number of excitons the sample can accommodate. 

As Frenkel excitons are constructed on highly localized electron-hole pairs, we could expect them to behave as 
elementary bosons. While their closure relation is indeed the one of elementary bosons - with a prefactor 1/N\ while 
the one for Wannier excitons has a prefactor (1/7V!) 2 - we find that they definitely behave as composite bosons through 
many other properties. In particular, the normalization factor of N Frenkel excitons (v\ BqBq 1 \v) = N\Fn differs 
from its elementary boson value N\ by a factor 

N 1 

Fn = (N S -N)W» (12 - 3) 

which, as for Wannier excitons, decreases exponentially with N-q. 

In addition to carrier exchanges, Frenkel excitons also interact through Coulomb processes. The interaction scatter- 
ing they produce, are derived, as for Wannier excitons, through a set of creation potentials. The part of the interaction 
scattering coming from direct Coulomb processes between sites depends on the momentum transfer as 

Zcoui($- + %%) = ^5qqiW q (12.4) 



Wq = £ e"«>- R [Vk (I \) + V R (g g) - V R (I ?) - V R (? I)] (12.5) 

R#0 

where Vr {^J, u J^j is the matrix element of the Coulomb potential between atomic sites at R, one electron going from 

the atomic state V\ to v' x while the other goes from v<i to v' 2 (see Eq. (2.6)). 

Besides this rather standard interaction scattering, we have identified a novel "transfer assisted exchange" scattering 
which has similarity with "photon assisted exchange" we found between polaritons. It comes from the indirect Coulomb 
processes responsible for the excitation transfer of Frenkel excitons, the coupling between two excitons being made 
through the Pauli scattering for carrier exchanges between these excitons. This "transfer assisted exchange" scattering 
depends on the "out" momenta (Q' 1; Q2) but not on the momentum transfer (Q^ — Qi) since the excitation transfer 
does not induce any momentum change. Its precise value reads 

Gran./ (g J) = (V Q1 + Vqj)A (§ J ) (12.6) 

where Vq is the Q dependent part of the exciton energy. 

Vq = E e ~ iQ R ^ (? o) (12-7) 

It is possible to show that, in this transfer assisted exchange, symmetry between "in" and "out" states is restored for 
energy conserving processes since we do have 

t (Q' 2 Ch\ _ p (Q2 Q' 2 \ 

Straus f QjJ Straus f \^q 1 J 



= {E Q[ + E Q , E Ql E Q2 ) X (J? J) (12.8) 
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Due to this "transfer assisted exchange" scattering, here also it is not possible to construct an effective bosonic 
Hamiltonian for Frenkel excitons which is hermitian, as physically required, and which produces the same matrix 
element as the one of the exact fermionic Hamiltonian. 

Frenkel excitons turn out to be quite tricky bosons: they appear as elementary in some cases but composite in many 
others. This is why it seems hard to guess physical results like the scattering rate of Frenkel excitons, the ground 
state energy of N Frenkel excitons or the nonlinear susceptibilities of materials having such excitons, without a precise 
calculation of these quantities using the procedure developed in this paper. These problems will be addressed in a 
near future. 
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Figure captions 

Fig 1. (a) Electron-hole pair eigenstate aJ^&J^ \v) of H e h defined in Eq. (2.3). 

(b) Electron-hole pair eigenstate a^b^ \v) of H pair defined in Eq. (2.4). 

(c) Indirect Coulomb potential V tra ns as defined in Eq. (2.10) which allows to transfer the excitation from site n-i 
to site n\. 

(d, c, h) Direct Coulomb potentials V e h, V ee and Vhh as defined in Eqs. (2.26-28). 

Fig 2. (a) Shiva diagram for the Pauli scattering A ^p, p^j of two "in" Frenkel excitons (P, Q) toward two "out" 

excitons (P', Q') as defined in Eqs. (3.11) and (4.2, 4.4). 

(b) Shiva diagram for carrier exchange of three Frenkel excitons defined in Eq. (4.5) and its decomposition in 2 x 2 
Pauli scatterings as given in Eq. (4.6). 

(c) Same as (b) for four excitons. 
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(d) Two consecutive carrier exchanges should reduce to an identity as seen from Fig. 2(e), while two consecutive 
Pauli scatterings do not reduce to this identity, as seen from Eq. (4.7). 

Fig 3. Shiva diagram for the scalar product of two Frenkel exciton states, as given in Eq. (5.2). 
Fig 4. Shiva diagram for the scalar product of three Frenkel exciton states as given in Eq. (5.5). 

Fig 5. (a) Part of the interaction scattering £ cou ; (^f associated to direct Coulomb processes between electrons 

and holes of different sites defined in Eq. (8.13). 

(b) Part of the interaction scattering £ trtms (^f associated to the indirect Coulomb processes insuring the 

excitation transfer defined in Eq. (8.14). This conceptually new scattering can be seen as a "transfer assisted 
exchange", similar to the "photon assisted exchange" scattering we have identified between two polaritons. 



